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ABSTRACT 

The Wide-field Infrared Survey Explorer has surveyed the entire sky at four in- 
frared wavelengths with greatly improved sensitivity and spatial resolution com- 
pared to its predecessors, the Infrared Astronomical Satellite and the Cosmic 
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Background Explorer. NASA's Planetary Science Division has funded an en- 
hancement to the WISE data processing system called "NEOWISE" that allows 
detection and archiving of moving objects found in the WISE data. NEOWISE 
has mined the WISE images for a wide array of small bodies in our Solar System, 
including Near-Earth Objects (NEOs), Main Belt asteroids, comets, Trojans, and 
Centaurs. By the end of survey operations in February 2011, NEOWISE iden- 
tified over 157,000 asteroids, including more than 500 NEOs and ~120 comets. 
The NEOWISE dataset will enable a panoply of new scientific investigations. 



1. Introduction 



A full understanding of the history of our Solar System requires study of its most 
primitive and dynamic members, the asteroids and comets. These different classes of small 
bodies are related; recent models of planetary migration, small body coUisional history, and 
the role that primitive bodies have played in shaping our planet's evolution predict that these 
signatures are imprinted on the numbers, distribution, transport and physical properties of 
these objects (Gomes et al. 2005 Tsiganis et al.||2005| [MorbideUi et al.||2005D . Only through 
detailed physical studies of all the major populations of our Solar System's small bodies can 
their origins, evolution and fate be learned. 

To date, the most prolific surveys to discover and characterize asteroids and comets have 
been conducted at visible wavelengths (e.g., the Catalina Sky Survey, the Lincoln Near-Earth 
Asteroid Research Program, the Near Earth Asteroid Tracking Program, and Spacewatch - 
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missions to primitive bodies is necessarily small (the Near-Earth Asteroid Rendevous mission, 
Giotto, Galileo, Deep Space 1, Deep Impact, Hyahusa, Dawn, Stardust, Rosetta Cheng|[2002 



Oberst et aI][200H |Reinhard|pM| [Thomas et"al]p^94l |A'Hearn et~ar]|2005| [Fujiwara et 



al. 2006 Russel et al. 2007 Brownlee et al. 2003; Barucci et al. 2005). Radar observations 



(Ostro et al. 2007) are a powerful means of obtaining detailed physical characterization 
remotely, but the objects must pass quite close to Earth. Hence, although approximately 
7600 near-Earth objects (NEOs), hundreds of thousands of Main Belt asteroids (MBAs) and 
~3000 comets have been discovered to date, detailed physical characterizations have been 



made for only a small fraction of these objects (e.g., Binzel et al. 2004; Wolters et al. 2008 



Harris et al. 2009; Bhattacharya et al. 2010). 



3 



1.1. Infrared Properties of Primitive Bodies 



Thermal infrared observations offer a powerful means of characterizing physical proper- 
ties of primitive bodies, and all-sky infrared surveys offer an opportunity to obtain these data 
for large numbers of objects in a relatively short time. Once an orbit has been calculated for 
an asteroid, its distance from the Earth at any given time is known, allowing the measured 
flux to be related directly to the flux emitted at the body's surface. Likewise, the known dis- 
tance to the Sun allows derivation of the incident solar flux. The physical surface properties 
of the asteroid can then be modeled to reproduce the measured emitted flux as a function 



of the incident radiation (Harris & Lagerros 2002 Lebofsky & Spencer 1989). Diameters 



derived from infrared measurements have smaller uncertainties than diameters derived from 
visible-light observations. Thermal flux emitted by an asteroid has a weaker dependency on 
a body's unknown albedo value (p^) than reflected sunlight. With both visible and IR data, 
the degeneracy between diameter and albedo can be broken, and both quantities can be 
independently fit. Current estimates for the range of albedos for near-Earth objects (NEOs) 
result in a factor of five uncertainty in diameter for a given reflected optical flux ( [Stuart fc 



Binzel 2004), since optical flux is linearly proportional to py. Yet the variation of an aster- 
oid's thermal flux is proportional to (1 — qpvt^^ , where q is the phase integral (Lebofsky & 



Spencer 1989). Strong constraints on diameter can still be obtained even for objects which 



have thermal fluxes but no visible measurements. For example, for a 0.13 km NEO with a 
well-known orbit observed to Sex at 12 /xm at a distance of 0.5 AU from Earth and 1.1 AU 



from the Sun, the diameter can be determined to within 10% (Wright et al. 2010 Mainzer 



et al. 2011) 



Because thermal flux is more weakly dependent on albedo than visible flux, infrared 
surveys are less subject to the bias against low albedo objects than visible light surveys. 



Although considerable work has been done to debias visible surveys (e.g., Stuart & Binzel 



2004 Jedicke & Metcalfe 1998; Harris 2008), IR surveys offer a complementary approach to 



determining the size and albedo distributions of small bodies in the Solar System. 

The Infrared Astronomical Satellite observed ~1800 asteroids, most in the Main Belt 
{IRAS; Tedesco et al.|1988 Matson, D.|1986 ). Ground-based thermal infrared measurements 
of asteroids have been made of ~100 objects ( Delbo et al.||2003 ). The Spitzer Space Telescope 
has undertaken a number of programs to observe Solar System objects, including 700 NEOs 



(Trilling et al. 2010). However, almost all objects observed by IRAS, Spitzer, and ground- 



based infrared telescopes have been optically selected, so the bias against low albedo objects 
is still present, complicating interpretation of the size and albedo distributions derived from 
these data. A large-area infrared survey capable of independent detection and discovery of 
a statistically significant number of objects is needed. With the Wide-field Infrared Survey 
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Explorer (Wright et al. 2010), such a resource is now available; the WISE bandpasses cover 
the thermal peak of most Solar System objects out to ~9 AU. The WISE data processing 
pipeline has been augmented to allow identification and discovery of new moving objects 
independently of previous visible asteroid surveys; these modifications are collectively known 
as "NEOWISE". Since it is capable of independent discovery of new objects, NEOWISE 
is unique in its ability to carry out a survey that is much less dependent upon albedo 
than other surveys. The purpose of this paper is to introduce the community to the rich 
NEOWISE dataset, including the pipeline that was used to detect moving objects, and to 
give an overview of the science enabled by NEOWISE. More detailed results will follow in 
separate papers. 



1.2. The Wide-field Infrared Survey Explorer 

WISE is a NASA Medium-class Explorer mission designed to survey the entire sky 
in four infrared wavelengths, 3.4, 4.6, 12 and 22 /im (denoted Wl, W2, W3, and Wi 
respectively) ( [Wright et aL]|2010[ |Liu et al.|[2008j [Mainzer et aL]|2005D . WISE is producing a 



multi-epoch image atlas and source catalog that will serve as an important legacy for future 
research. The mission's two primary science goals are to find the most luminous galaxies in 
the entire universe and to find the closest and coolest stars; a secondary goal is to study the 
nature of Earth's nearest neighbors, the asteroids and comets. 

WISE consists of a 40 cm telescope in a Sun-synchronous low-Earth orbit. The telescope 
scans the sky continuously, pointing near the zenith while a scan mirror freezes the sky on 
the focal planes for approximately 11 seconds. While the sky is frozen on the focal plane, an 
8.8 sec exposure is taken of the same 47' x 47' field of view in the four WISE bands. The 
scan plane precesses by ~l°/day. Each image overlaps by 10% in the in-scan direction, with 
90% overlap in the cross-scan direction. This orbit allows WISE to cover nearly every part 
of the sky a minimum of eight times, with more coverage at the ecliptic poles. The survey 
began on 14 January 2010, and full coverage of the sidereal sky was achieved in mid- July, 
2010. The telescope and detectors were maintained at their operating temperatures through 
the use of a dual-stage solid hydrogen cryostat. The secondary cryogen tank depleted on 5 
August 2010, and the primary tank was exhausted on 30 September 2010. After that time, 
the first two channels {Wl and W2) continued to operate, and a four month survey (known 
as the NEOWISE Post-Cryogenic Mission) using these two channels only and the same 
survey strategy was undertaken. The NEOWISE Post- Cryogenic Mission was completed on 
1 February 2011. 



Astrometric errors are less than 0.5 arcsec with respect to 2MASS for objects with high 
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signal-to-noise measurements. The preliminary estimated SNR=5 point source sensitivity 
on the ecliptic is 0.08, 0.1, 0.85 and 5.5 mJy in the four bands (assuming eight individual 



exposures per band; Wright et al. 2010). Sensitivity improves away from the ecliptic due to 



denser coverage and lower zodiacal background. For an asteroid with a bolometric albedo 



of 0.05, the W3 sensitivity of 0.85 mJy corresponds to an visible magnitude of 23 (Wright 



et al. 2010). However, as described below, the search for new minor planets was carried out 



using individual exposures. 

The WISE survey cadence yielded excellent repeated coverage of moving objects in the 
Solar System, with an average of 10 detections of a typical asteroid or comet spaced over 
~36 hours. This is a slightly higher number of detections than the corresponding value of 
eight for an inertially fixed object on the ecliptic plane because asteroids tend to move in 
the same direction as the WISE scan plane. 

As of February 2011, WISE detected more than 157,000 moving Solar System objects. 
The dataset contains at least 584 NEOs, ~2000 Jupiter Trojans, 18 Centaurs and scattered- 
disk objects, and ~120 comets. WISE has also detected a handful of irregular satellites of 
Jupiter and Saturn, including lapetus and Phoebe, as well as cometary dust trails and dust 
bands in the zodiacal cloud. 



2. NEOWISE: An Enhancement to the Baseline WISE Mission 

Data processing in the baseline WISE mission does not include any provision for iden- 
tifying previously unknown Solar System objects. Furthermore, only the coadded Atlas 
images that combine many exposures covering each point on the sky were planned to be 
pubhcly released. Outlier rejection algorithms designed to reduce susceptibility to transient 
phenomena such as cosmic rays suppress most moving objects from the coadded images. 
To facilitate Solar System science with WISE, two enhancements were made to the baseline 
pipeline; these tasks are collectively called NEOWISE ("Near Earth Object + WISE"). 



2.1. NEOWISE Task 1: Creation of the archive of single-exposure images plus 

searchable interface 

NEOWISE Task 1 allows researchers to query the WISE data archive for pre-discovery 
astrometry and physical data for Solar System objects that are discovered after the final 
data processing is complete. It also enables users to download the individual single-epoch 
images of these objects. 
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Recovery of small bodies in the WISE data archive requires long-term access to the 
time-tagged, cahbrated images and the database of extracted sources from the individual 
(or single-epoch) WISE exposures. These are referred to as the Level lb image and source 
products. The Level lb images are FITS format images in each of the four bands that have 
instrumental signatures removed (e.g. linearized, droop-corrected, dark-subtracted, flat- 
fielded, illumination corrected, etc) and that have the reconstructed WCS and photometric 
zero point information in their headers. The Level lb extracted source database contains 
the reconstructed equatorial positions, calibrated brightness in four bands and associated 
uncertainties, and assorted quality and characterization flagging for sources detected on the 
Level lb images. Latent images and other instrumental artifacts are flagged in the source 
database. These products are produced by the Scan/Frame Pipeline component of the WISE 
Science Data System (WSDS). 

The Atlas, Catalog and assorted metadata products are archived for the long-term and 
served to the research community and public via the online and machine- friendly interfaces of 
the NASA Infrared Science Archive (IRSA). An enhanced interface to the IRSA catalog and 
image query services that allows searches based on Solar System object name (via Horizons 
name resolution) or orbital elements has been developed as part of NEOWISE Task 1. 



2.2. NEOWISE Task 2: Rapid Identification and Posting of Previously 
Unknown Solar System Object Candidates 

NEOWISE Task 2 consisted of an augmentation to the ground data processing system 
that exploits the WISE observing cadence to identify moving objects by comparing source 
extractions made from imaging of the same region of sky on successive orbits. New moving 
object candidates were reported to the Minor Planet Center (MPC) within ten days of 
the midpoint of their observations on board the flight system. The Minor Planet Center 
serves as the central clearinghouse for NEOWISE astrometry and associates it with all other 
observations. 



2.3. WMOPS: The WISE Moving Object Processing Software 



Candidate moving object tracklets (sets of time-tagged positions) were identified using 
the WISE Moving Object Processing System (WMOPS). WMOPS was developed using 



elements of the PanSTARRS Moving Object Processing System (MOPS; Denneau et al. 



2007 Jedicke et al. 2009) that linked detections of candidate moving objects from multiple 
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exposures covering the same regions of the sky. WISE-specific elements of WMOPS included 
optimizing detection for the WISE survey cadence, adapting to WISE detector and source 
extractor characteristics, and specialized quality assurance filters and tools to allow vetting 
of moving objects prior to delivery to the MPC 

WMOPS was integrated into the general WISE data processing pipeline system and 
followed scan/frame processing. The scan/frame pipeline performed instrumental image 
calibration (linearization, dark-subtraction, flat-fielding, sky subtraction, etc.) on individual 
WISE images from one scan (where a scan contains ~260 four-band exposures taken during 
1/2 of a WISE orbit), detected and measured position and fluxes for sources in the calibrated 
images, and flagged sources that may be spurious detections of image artifacts produced by 
bright stars. Scan/frame processing produced photometrically and astrometrically calibrated 
single-exposure (Lib) images and extracted source lists. 

WMOPS processing steps are illustrated in Figure [Tj Input was the source lists from 
~90-120 scans (spanning approximately 3-4 days). The input source lists were first filtered 
to exclude low reliability extractions such as flagged artifacts, very low signal-to-noise ratio 
sources, and extractions with pathologically large values of the PSF-fitting chi-squared val- 
ues. Stationary object rejection (SOR) was then performed by positionally associating the 
filtered source lists from overlapping images and identifying pairs of detections that do not 
move by more than 2 arcseconds between adjacent exposures, and eliminating these detec- 
tions from further processing. After SOR, the remaining moving object candidate detections 
are processed by the findtracklets routine that creates pairs of detections within a region of 
sky and correlates each pair's velocity and direction of motion, along with their sky location. 
The pair lists are stored in a database and then fed through the collapsetracklets routine 
that creates chains of pair-correlations by matching pairs according to velocity and position 
within a specified uncertainty threshold at each point. Candidate tracklets were required to 
have a minimum of five independent detections to be considered reliable. 

WMOPS efficiently detected all moving objects whether they were previously known or 
unknown. Approximately 80% of the tracklets were previously known Solar System objects 
and were automatically accepted. The remaining 20% of tracklets were validated by visual 
inspection by NEOWISE scientists. Figure |2] illustrates the web-based quality assurance 
tool that enabled rapid assessment of a tracklet's reliability. Tracklets for known objects 
and reliable new detections were reported to the MPC. 

WMOPS was run approximately every four days to ensure that the accumulated po- 
sition uncertainties for the fastest moving objects such as NEOs remained small enough 
to enable efficient recovery by ground-based observers. Tracklets from each WMOPS run 
were sent to the MPC for publication and linkage to other observations. Objects that the 
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MPC determined to have a high probabihty of being NEOs or comets were posted to the 
NEO confirmation page ( jhttp: / /www.min orplanetcenter.org/iau/NEQ /ToConfirm.html ) to 
be followed up by the world-wide network of professional and amateur astronomers. Ground- 
based follow-up was essential for NEOs in particular because the WISE observational arcs, 
while generally long enough to assign a designation to new objects, do not span a sufficient 
timebase to allow recovery at subsequent apparitions. Ground-based follow-up within two 
to four weeks allows well-determined orbits to be computed. Objects that are confirmed to 
be NEOs or comets are published by the MPC via the Minor Planet Electronic Circulars 
(MPECs) or lAU Circulars. The astrometry and computed orbits of slower moving objects 
were published on a monthly basis in the MPC's Minor Planet Orbit Supplement. 



2.3.1. Range of Motion Limits 

The minimum tracklet length validation requirement and the minimum sky-plane pro- 
jected velocity imposed by the SOR defines the range of motion to which WMOPS is sen- 
sitive. Parallax motions below 2 arcsec over 90 minutes (approximate time between WISE 
coverages near the ecliptic) correspond to a distance of roughly 70 AU. The limits of five 
detections imposes an upper limit in that objects may outpace the survey cadence before 
the minimum number of detections can be obtained. Motions greater than the upper lim- 
its perpendicular to the scan direction are roughly 3.25°/day, considering both an irregular 
cadence pattern as well as an evenly progressing survey. The maximum speed of an object 
detected by WMOPS was 3.3°/day, but the component perpendicular to the scan direction 
was 3.227day. 



3. Results 

Figures 4a and 4b show the distribution of orbital elements of all the Solar System 
objects identified by NEOWISE to date, including NEOs, MBAs, comets, Trojans and Cen- 
taurs. By February 2011, more than 252,000 tracklets were reported to the MPC during the 
cryogenic and post-cryogenic missions, representing observations of ~ 157,000 unique minor 
planets. 
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3.1. Near-Earth Objects 

NEOs are thought to be a transient population, persisting in their orbits for only a 
few million years before being ejected into the outer Solar System, crashing into the Sun, 



or colliding with other bodies (Gladman et al. 2000). Therefore, they must be continually 



resupplied from more permanent source reservoirs such as the Main Belt and comets. Recent 
models have improved our understanding of the transport mechanisms and sources of NEOs 



(e.g., Bottke et al. 2002) based on observations of their size frequency distribution (SFD). 
However, current NEO SFD models rely on the absolute magnitude {H) through an assumed 
albedo to derive size. By observing and discovering a sample of NEOs at thermal infrared 
wavelengths, we can obtain an improved computation of the SFD because 1) IR-derived 
diameters are intrinsically more accurate; 2) we can obtain albedos and further constrain 
diameters for all individual objects that have measured optical magnitudes; and 3) an IR 
survey is less biased against low albedo objects. The improved SFD offered by NEOWISE 
will allow us to test NEO population models, to study the coUisional histories of NEO and 
MBA populations, as well as to examine how SFD and albedo vary with asteroid collisional 
families. 

Our understanding of the hazard posed by NEOs also depends on knowledge of their 
SFD, and efforts to plan mitigation strategies depend on the physical properties of NEOs. 
To the extent that albedo can be used as a proxy for asteroid composition, it will be possible 
to derive an approximate density, mass, and estimated impact energy for previously known 
and new objects observed by WISE. 

As of February 2011, NEOWISE has identified more than 584 NEOs, including discov- 
ery of 135 previously unknown NEOs. Figure [5] shows the preliminary distribution of H 
magnitudes for all NEOs known to have been observed by NEOWISE during the cryogenic 
mission. The distribution peaks at /7~18, corresponding to diameters between ~470 and 
1500 m for albedos between 0.05—0.5 {H values for previously known objects were taken 
from the MFC observation file; H values for newly discovered objects were obtained by visi- 
ble follow-up observations from a worldwide network of amateur and professional observers). 
The survey cadence includes many repeated coverages at the ecliptic poles, regions which 
most ground-based surveys preferentially avoid. NEOWISE is thus sensitive to objects at 
extreme declinations, which are more likely to have high orbital inclinations or have orbits 
that bring them very close to the Earth (Figure |6]). Furthermore, unlike ground-based tele- 
scopes, the survey is unimpeded by daytime, weather, or variable seeing, so it is easier to 
debias (c.f. Jedicke fc Metcalfe||1998 ). However, of the 135 new NEOs discovered, ~15 NEOs 



were designated by the MFC but received no optical follow-up either due to their optical 
faintness or unavailability of suitable ground-based facilities. Due to the ~36 hour arc length 
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typical of most NEOWISE detections, ground-based follow-up was necessary in most cases 
to obtain a secure orbit for the newly discovered NEOs. An additional ~15 NEO candidates 
were placed on the MFC's NEO Confirmation Page but received neither designations nor 
follow-up. A radiometric diameter can be determined for these objects (subject to the un- 
certainties in Sun-object-observer distances due to orbital uncertainties), but albedo cannot 
be computed without optical fluxes. It is hoped that these objects will be recovered either 
from archival data or future surveys. Without this, these objects will result in an additional 
uncertainty in the derived NEO albedo distribution. 



3.2. Main Belt Asteroids 



NEOWISE has observed over 154,000 MBAs as of February 2011, including ~33,000 
new discoveries. This represents an advance of nearly two orders of magnitude over the 



IRAS Minor Planet Survey (Tedesco et al. 2002), allowing a correspondingly more detailed 



view of the size and albedo distributions throughout the Main Belt. Figure [7] shows that the 
preliminary distribution of H magnitudes for all MBAs observed by WISE by the end of the 
cryogenic mission peaks at ~15, corresponding to diameters between 1.9—5.9 km for albedos 
between 0.05—0.5. H values were taken from the MFC observation files for previously known 
objects; however, the H values for the newly discovered objects are approximated from the 
objects' thermal fluxes since no dedicated visible follow-up for these objects was carried 
out. It is expected that linkages will eventually be made between visible survey data and 
WISE observations, resulting in increasing numbers of NEOWISE-discovered objects with 
well-measured H values and refined orbits. Figure [8] shows the preliminary H magnitude 
distribution as a function of semimajor axis within the Main Belt. With thermal infrared 
fluxes available for all of these Main Belt asteroids, it is now possible to derive diameters and 
(where optical fluxes are available) albedos. This will facilitate study of the distribution of 



sizes and albedos within dynamical families (Hirayama 1918) as has been done with multi- 



band optical photometry (Parker et al. 2008 Ivezic et al. 2002) 



The NEOWISE dataset contains observations of asteroids that have been previously 



associated with one of the 54 dynamical families defined by Nesvorny et al. (2006a). Figure 



|9] shows the number of objects observed by NEOWISE for each family; for some families, 
NEOWISE has observed thousands of known members. The combination of WISE infrared 
fluxes with multi-band optical photometry should allow new family members to be identified 



on the basis of their albedos and optical colors as well as their dynamical properties. Parker et 



al. (2008) have shown that it is possible to use Sloan Digital Sky Survey (SDSS) colors to find 



two to three times more family members that make up "halos" surrounding existing families. 
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These are family members that have moved far enough into the background population that 
they can no longer identified by hierarchical clustering methods, yet they have colors identical 
to the cores of the families (e.g., Vesta and Eos). In particular, this technique may facilitate 
the association of the smallest, least massive objects with various families. These objects 
are of interest because their small sizes make them more susceptible to the non-gravitational 
forces that cause the most spread in the orbital elements which are used to compute family 
ages (c.f. Nesvorny et al. 2006a). By incorporating size and albedo distributions derived 
from infrared fluxes, as well as by finding additional family members using their albedos and 
colors, the NEOWISE dataset should facilitate an improved understanding of the ages and 
timing of the breakups that created the present-day Main Belt families. 



3.3. Comets 



As reservoirs of some of the most pristine primordial material in the Solar System, 
comets serve as Rosetta stones for understanding the transport and history of volatiles. 
Comets sample the most distant regions of the Solar System: the Kuiper Belt just outside of 
Neptune's orbit, and the Oort Cloud, which may extend halfway to the nearest star. While 
we cannot directly observe the earliest stages of comet formation in our own Solar System, 
the size frequency distribution and orbital distribution of present-day comets preserve a 
record of the outer solar nebula's mass distribution that can be used to model the coUisional 



and dynamical evolution of comet source regions over the last 4.5 Gy (c.f. Tsiganis et al. 



2005 Levison et al. 2009). The cometary size and orbital element distributions allow a test 



of models of the origins of Halley family comets (HFCs) and long period comets (LPCs) 
as well as Jupiter family comets (JFCs). Low-inclination JFCs are thought to originate in 
the Kuiper belt and are expected to have smaller sizes than the HFCs and LPCs due to 



coUisional processes and mass loss caused by their numerous passages by the Sun (Weissman 



& Levison 1997). The HFCs and LPCs, believed to originate in the Oort cloud, should 



have retained their original size distribution (Stern & Weissman 2001) and should be larger 



(Meech et al. 2004). Further, the distribution of perihelia can set constraints on the number 



and mass of objects in the Oort cloud (Francis 2005). 



WISE has observed ~120 comets as of February 2011, of which 20 are new discoveries 
(including three previously known asteroids on which WISE observed activity; these obejcts 
were subsequently redesignated as comets). The NEOWISE survey of comets offers a dataset 
that is less biased against low-albedo objects than previous optical surveys, and as an all- 
sky survey, it has covered high declination regions as well as the ecliptic plane with a well- 
determined set of survey pointings. Detections of previously known comets were identified by 
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WMOPS, and new discoveries were made by inspecting the ~20% of tracklets that were not 
associated with known objects to look for evidence of extended emission. New discoveries 
suspected of being comets were posted to the MFC's NEO Confirmation Page. All of the 
20 new discoveries received visible follow-up, allowing orbits to be refined beyond the initial 
orbits computed solely from NEOWISE observations. As with the other major populations 
detected by WISE, infrared fluxes can be used to compute nucleus size when cometary 
activity does not significantly obscure it. While the Spitzer Space Telescope has undertaken 



a survey of 100 comets (Fernandez et al. 2008), they are all JFCs that have been selected 



by optical surveys. The NEOWISE survey includes both long and short period comets, 
as well as independent discoveries, allowing constraints to be set on the ratio of long to 
short period comets. Of the ~120 comets observed by the end of the cryogenic portion 



of the mission, 35% are LPCs. Figures lOab show the distribution of orbital elements of 



comets detected by NEOWISE during the cryogenic mission. Many of the comets have high 
inclination orbits; however, the survey cadence resulted in increased coverage at the ecliptic 
poles, leading to a bias in favor of detecting comets at high inclinations (or closer to Earth). 
Of the 20 comets discovered by NEOWISE as of February 2011, approximately one-third 
have perihelia greater than 4 AU, and the remaining two-thirds are LPCs, four of which have 
nearly parabolic orbits. 

For objects where significant cometary activity is observed, it is possible to compute 
particle size distributions, allowing the strength of non-gravitational forces to be determined 
as a function of particle size. For comets with sufficiently high fluxes at bands Wl (3.4 yum) 
and W2 (4.6 /im), gas ratios of CO and CO2 can be computed. 

NEOWISE has also observed four objects that have been targeted by spacecraft for 
in situ visits, 19P/Borrelly {Deep Space 1), 9P/Tempel 1 {Deep Impact), 103P/Hartley 



(EPOXI; Bauer et al. 2010) and 67P/Churyumov-Gerasimenko, which the Rosetta mis- 



sion will encounter in 2014. Figures llabc show the WISE views of 67P/Churyumov- 



Gerasimenko at two apparitions. A preliminary flt of a blackbody function to the flrst 
epoch of observations on 19 January 2010 yields an effective temperature of 179 K using an 
aperture radius of 22 arcsec. A similar flt to the second epoch observations on 4 July 2010 
yields 155 K (Figure 11c). These fltted temperatures exceed ambient blackbody tempera- 
tures by about 20 K. It is not clear that that this is due to the dust grain temperature, as 
there is a silicate emission band near the center of WISE band W3 at ~ 10-12 /xm. 
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3.4. Jupiter Trojans 



Trojans are populations of minor planets that are co-orbital with a planet. The known 
Jupiter Trojans (there are also a handful of objects that share orbits with Mars and Neptune) 
are separated into two clouds of objects that librate around the L4 and L5 Lagrangian points 
in Jupiter's orbit. These two points are found 60° of heliocentric ecliptic longitude ahead 
(also called the leading cloud) and behind (the trailing cloud) the planet, and both clouds 



possess enough dynamical stability to survive over the age of the Solar System (Marzari & 



SchoU 1998). Reviews of the Trojan populations can be found in Jewitt et al. (2004) and 



separately in Dotto et al. (2008 ). The Trojans are an important population for understanding 



the origin and evolution of Jupiter, and they also provide a potential window into the solar 
nebula at the distances of the giant planets, since they are the only major population that 
might have survived the dynamical clearing of minor planets in this region during the early 
formation of the Solar System. But how and when the Trojans were trapped in Jupiter's 
1:1 mean-motion resonance remains unknown. Whether the capture happened during a very 
early epoch while the planet was forming or much later during the dynamical clearing phase 



is still under discussion (Morbidelli et al. 2005 Marzari & SchoU 2007). 



There are currently about 4700 known Jupiter Trojans, and NEOWISE has reported 
observations of ~1600 of these, increasing the number of thermal emission measurements of 



this population by more than an order of magnitude compared to earlier surveys ( Cruikshank 



et al.|1977t[Tedesco et al.|2002| [Fernandez et al.|2003|[2009| . In addition, NEOWISE detected 
an additional 400 objects that the MFC has designated with Trojan-like orbits, but these 
have arcs that are too short for their identification to be secure. WISE observes the sky in 
scans of constant ecliptic longitude, and as the survey began, the first scans were already 
in the middle of the leading Trojan cloud. The scans moved toward the tail of the cloud, 
resulting in about 20-30% of the leading cloud closest to Jupiter remaining unobserved during 
the cryogenic survey (although it was observed about seven months later with bands Wl 
and W2 during the post-cryogenic mission). The survey exited the leading cloud at the end 
of February, 2010 and almost immediately entered the trailing cloud on the opposite side of 
the Solar System (Figure [3]). Observations of the trailing cloud were completed in late May, 
2010. The tail end of the trailing cloud was close to the galactic center when it was observed 
by NEOWISE, leading to noticeable loss of detection efficiency for less than 10% of the 
trailing cloud. While the two clouds were not observed with perfect uniformity, NEOWISE 
has provided the most uniform survey of the Jupiter Trojan population to date. Of the 
~1600 objects that have been conclusively identified as Trojans, ~900 were observed the 
leading cloud and ~700 in the trailing cloud. This yields a non-debiased ratio of the number 
of objects in the L4 and L5 cloud of ~1.3. The 400 probable Trojans are distributed roughly 



equally among the two clouds. The H-distribution (Figure 12) shows that the survey's limit 
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is H~13 (or ~12 km for an assumed albedo of 0.07). 



The survey provides a unique opportunity to not only characterize the physical proper- 
ties of a large fraction of the Trojans, but it also provides a uniform survey of both clouds. 
The survey must be debiased in order to test the initial results that show an over-abundance 
of objects in the leading cloud compared to the trailing cloud. In addition, a proper vetting 
of the short arc sample using the on-sky velocity (similar to the method used on the SDSS 
data by Szabo et al. 2007) and thermal colors is needed to define the Trojan sample. As 



with the Main Belt asteroids, NEOWISE detections of new Trojan candidates may be linked 
to other observations in the future, allowing refinement of their orbits. 



3.5. Unusual Objects and Outer Solar System Objects 

WISE^s wavelengths are not long enough to detect large numbers of objects in the outer 
Solar System via their thermal emission. Nonetheless, several dozen objects were observed 
and discovered in this region with the NEOWISE pipeline. A number of unusual objects 
with high inclinations have been discovered by NEOWISE such as 2010 OAlOl, an outer 
Main Belt Asteroid with an inclination of 84°, and 2010 CR140, an object with a semi-major 
axis of 5.6 AU and an inclination of 74.7°. 2010 CR140 may represent an object in transition 
between the Jupiter Trojan and Centaur populations. In total, eleven objects were observed 



in the outer Solar System with inclinations larger than 80°. Figure 13 shows an a-e plot 
of WISE-observed objects, highlighting the outer Solar System. WISE has observed ~20 
Centaurs to date, including five new discoveries. 



3.6. Lightcurve Studies and Thermophysical Modeling 



By obtaining an average of ten observations spaced uniformly over ~36 hours for each 
Solar System object, the WISE cadence lends itself to the study of small body lightcurves. 
Ground-based observers are interrupted by daytime, clouds, and seeing variations; however, 
WISE is able to provide an extremely uniform, stable platform with regular observations 
spanning days or even, in a few cases, weeks. An example of this type of long-period 



lightcurve for a near-Earth asteroid is shown in Figure [14} This makes the WISE dataset 
unique in comparison with ground-based lightcurve surveys that are generally limited to 
studying rotational periods faster than ~1 day (Masiero et al. 2009 Pravec et al. 2008 



2006). The WISE observations offer a unique window into long-period lightcurve studies. 



and with suitable debiasing, the fraction of objects with rotational periods longer than a 
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day may be determined. Figure 15 shows the distribution of the number of observations for 
objects detected by WMOPS; on average, a typical Solar System object was detected every 
three hours by NEOWISE. 

Binary asteroids have been found via photometric observations, direct imaging (ground- 
based and space-based), adaptive optics, and radar. Binary systems provide a useful labora- 
tory because they allow masses to be determined; when combined with infrared diameters, 
density and porosity may be computed. It is possible to search for binary systems in the 
WISE dataset by searching for objects with lightcurve amplitude variations ^1.5 magni- 



tudes; such objects are the most likely to be binary systems ( Chandrasekhar 1969 Holsap- 



ple 2001). An example of a Main Belt asteroid displaying the large-amplitude lightcurve 



indicative of a candidate binary system is shown in Figure 16 With a period of more than 
2.5 days and an amplitude of ~1.5 magnitudes, this object, (49137) 1998 SC35, will require 
a concerted follow-up observing campaign to confirm its binarity. Preliminary fits to the 
individual WISE data points using a thermophysical model assuming a single spherical ob- 



ject and the Near-Earth Asteroid Thermal Model (NEATM; Harris 1998) yield diameters 
ranging from 3.8 to 7.5 km; the wide range of diameters is due to rotational effects and is 
much larger than the measurement uncertainties. 

As with visible data, shape and rotation state can be derived from WISE observations 



c.f. Delbo et al. 2003 Wright et al. 2007). If sufficient data are available at a range of 



different phase angles, it is possible to compute the pole orientation and rotation state of an 



object in addition to its shape model through lightcurve inversion (Kaasalainen et al. 2002 



Miiller et al. 



2005 



Durech et al. 2010). However, unlike models computed with visible data 



alone, the infrared fluxes allow the absolute size of the object to be determined as well. For 
some WISE-ohseived objects, it may be possible to derive an albedo map across an object's 
surface. 



3.7. Comet Trails and Dust Bands 



WISE thermal observations are sensitive to extended structures related to the principal 
mechanisms of dust production in the inner Solar System: comet dust trails and zodiacal dust 
bands. Dust trails are narrow (arcminute) structures extending degrees to tens of degrees in 
the sky, and they consist of low-velocity millimeter to centimeter diameter particle emissions 
from short-period comets ( Sykes et al.||1986a Sykes &: Walker||1992 ). Spitzer Space Telescope 



observations have shown that 80% of the short-period comets surveyed have trails (Reach 



et al. 2007). A cratering or disruptive coUisional event between two asteroids would also 



result in the formation of a dust trail, though the ejection velocities (compared to trail- 
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forming comet emissions) are expected to be higher, resulting in a broader trail (e.g. Sykes 



1988 Nesvorny et al. 2006b). Over time, gravitational perturbations acting on the trail 



particles results in the differential precession of their nodes, resulting in the formation of 
a torus of collisional material extending around the Sun. Some of these bands have been 
associated with the breakups of asteroids in the Main Belt that lead to the formation of 
young asteroid families ( Nesvorny et al.|[2003 , 2006b). Identifying the dust bands associated 
with young families will serve to confirm their age and will allow quantitative studies of 
the dispersal time for the bands (Espy 2009). WISE has already detected a number of 



unn (Figure 


17 


), which 


Colangeli et al. 


1998 


)• 



WISE image is approximately a hundred times more sensitive than IRAS at 22 /im and has 
significantly better spatial resolution. Future work should allow discovery of more trails that 
can potentially be associated with either asteroid families or terrestrial meteor streams. 



4. Conclusions 



With NEOWISE, we have detected over 157,000 minor planets during a year-long survey, 
including discovery of ~33,000 new minor planets. We have shown that the WMOPS system 
is capable of rapid identification of moving objects over a wide range of apparent velocities 
when provided with appropriately calibrated input images and source detections. The unique 
nature of the NEOWISE dataset will facilitate a wide range of small body studies from fields 
as diverse as population studies, infrared lightcurve shape modeling, thermophysical studies, 
comets, and dust bands and comet trails. With close to 100 times more objects observed than 
its predecessor, IRAS, NEOWISE will leave an even greater legacy for future researchers. 
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Fig. 1. — The flowchart summarizing WMOPS. 




advance to: 
O From Page 
©Next Object 

Fig. 2. — A sample of the quality assurance pages used to validate a tracklet for a previously 
unknown moving object. The top two rows show cutouts of the images in bands W3 and 
WA, respectively; the following rows show W3 images of the same inertial position on the sky 
before and after the exposures containing the target. These rows allow verification that the 
source observed is not inertially fixed. The final two rows show the exposures immediately 
preceding the top row; this allows verification that the transient source observed in the top 
frame is indeed a moving astronomical object and not an artifact fixed in pixel space but 
not in time. 




Fig. 3. — WISE always surveys near 90° solar elongation; this figure shows a top-down 
view of the objects detected by NEOWISE as of February 2011 (distances are given in AU). 
The outermost circle represents Jupiter's orbit; the interior circles represent the terrestrial 
planets. Previously known NEOs are shown as green circles; new NEOs discovered by 
NEOWISE are shown as red circles; previously known comets observed by WISE are shown 
as cyan squares, and comets discovered by NEOWISE are shown as yellow squares. All other 
objects are shown as black points. The drop in density of objects observed near (+2, +2) 
AU in the figure is due to the exhaustion of the secondary tank's cryogen on 5 August, 2010, 
resulting in the loss of band WA. The dashed red line indicates the exhaution of the primary 
tank and the start of the NEOWISE Post-Cryogenic Mission on 1 October 2010; the survey 
was completed on 1 February 2011. The drop in detections near (+2, -2) AU in the figure is 
due to the intersection of the galactic plane with the ecliptic plane; the higher backgrounds 
and confusion caused by galactic cirrus resulted in the identification of fewer sources. 
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Fig. 4a. — NEOWISE has detected objects throughout the Solar System. The orbital ele- 
ments of all objects that NEOWISE has detected in the inner Solar System with arcs longer 
than 30 days are shown as black dots; objects with arcs less than 30 days are shown in cyan. 
The quantization seen in the cyan points is due to the method by which the Minor Planet 
Center computes orbits for objects with short observational arcs. 
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Fig. 4b. — a-i plot of objects observed by NEOWISE in the inner Solar System; color scheme 
is the same as in Figure 4a. 
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Fig. 5. — The H magnitude distribution for the >500 NEOs detected by WISE as of Febru- 
ary, 2011. 
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Fig. 6. — WISE^s survey cadence results in increased coverage at the echptic poles, making it 
easier to detect NEOs at high inclinations. The preliminary, non-debiased NEO inclination 
distribution is shown here for all NEOs observed by WISE. The dataset must be debiased 
before conclusions about the NEO inclination distribution can be drawn. 
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Fig. 7. — The preliminary distribution of absolute magnitudes of ~150,000 Main Belt aster- 
oids detected by NEOWISE peaks at H ~15; this corresponds to diameters between 1.9—5.9 
km for albedos of 0.5 to 0.05. 
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Fig. 8. — The peak of the Main Belt asteroid H magnitude distributions shifts as a function 
of semimajor axis. The plot shows the distributions of four categories of Main Belt asteroids, 
as defined by Zellner et al. (1979): MBA I (2.06 < a < 2.5); MBA II (2.5 < a < 2.82); 
MBA III (2.82 < a < 3.27); and MBA IV (3.27 < a < 3.65). The solid line represents the 
histogram of MBA I, the dotted line represents MBA II, the dashed line is MBA III, and 
the dash-dot line is MBA IV. The histogram numbers of MBA IV have been multiplied by 
ten for clarity. 



Fig. 9. — Approximately 26,000 of the Main belt asteroids observed by NEOWISE during the 



cryogenic mission are objects that are known to belong to dynamical families Nesvorny et al. 



(2006a). In addition to these objects, it is highly likely that other previously unrecognized 
family members are distributed among the more than 150,000 other MB As detected by 
NEOWISE. These objects may be identified by their IR-derived albedos and optical colors. 
This will facilitate a wide range of studies aiming to improve the precision of calculations of 
parent body breakups. 
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Fig. lOab. — The preliminary, non-debiased distribution of orbital elements of NEOWISE- 
detected comets reveals many with high inclinations. 




Fig. llabc. — Rosetta target 67/P Churyumov-Gerasimenko. 11a, b: Three-color WISE 
images of the comet observed at two epochs: 19 January 2010 (left panel a) at a heliocentric 
distance R = 3.32 AU and an observer distance A=3.14 AU, and 29 June 2010 (right panel 
b) with R = 4.18 AU, A = 3.97 AU. In these images, the blue channel is band W2 (4.6 
/xm), green is band W3 (12 fim), and red is band W4 (22 fim). 11c: Blackbody fits to 67/P 
for the two epochs of WISE observations. Preliminary corrections to fluxes were made by 
interpolating between the blackbody color corrections provided in Wright et al. (2010). The 
blackbody fit for the first epoch of observations is shown as the solid black line, and the 
WISE data points are red diamonds; the blackbody fit for the second observational epoch is 
the dashed black line, and the data points are blue triangles. 67/P was undetected in bands 
Wl and W2 in both epochs. 
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Fig. 12.— The NEOWISE distribution of observed Trojans peaks at H ~13. Of the 4237 
Trojans known, WISE has observed ~1600 during the one year survey (sohd black hne), along 
with ~400 previously unknown Trojans for which orbits are presently poorly constrained. 
More Trojans were observed in the L4 cloud that leads Jupiter (red dashed line) than the L5 
trailing cloud (blue dash-dot line), despite the fact that the cryogenic portion of the survey 
ended before the survey of the L4 cloud was complete. 



i.or^ 



0.8 



t 0.6 
c 

0) 

B 0.4 

LU 



0.2 



0.0 



10 20 30 40 50 60 70 
Semi-major Axis (AU) 



Fig. 13. — The a-e distribution of objects observed by NEOWISE in the outer Solar System. 
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Fig. 14. — The 12 fim time series of the near-Earth asteroid 230118 (2001 DB3) extends 
over two days. The rotational modulation is readily visible due to WISE's unique observing 
cadence. 
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Fig. 15. — On average, WMOPS detected most moving Solar System objects ~10-12 times 
over 1.5 days; however, some NEOs have been observed more than a hundred times. 
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Fig. 16. — The 12 fim lightcurve of Main Belt Asteroid 49137 extends over 2.5 days and 
shows an amphtude variation large enough (~1.5 magnitudes) to indicate that the object 
may be a binary. The bottom of the lightcurve shows a flattened region that may indicate 
the signature of the eclipsing secondary. The object was detected by WMOPS beginning on 
Modified Juhan Date 55306.640. 
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Fig. 17. — The Jupiter-family Comet 65P/Gunn was observed by WISE on 24 April 2010. 
This single-frame observation shows 3.4 fim as blue, 4.6 /xm as green, 12 fim as orange, and 
22 fim as red. The nucleus and tail of the comet are clearly visible, along with the trail of 
debris that has been spread around the comet's orbit. 



